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^ ' Abstract 

: 

^ A global picture is drawn tying together most exact cosmological so- 

, lutions of gravitational theories in four or more spacetime dimensions. 

a^ ■ 
o : 

m . 1 Introduction 

O ■ 

Q , The search for exact solutions of gravitational theories in four spacetime di- 

mensions has often become an end in itself, a game in which one seeks to write 
}L{ '. down closed form expressions or at least quadratures which solve some particu- 

lar case of some set of equations, without worrying about interpretation, placing 
^ , them into context or drawing some significant conclusion about properties of the 

• original equations. Now this game, increasingly more difficult to play in four 

rS I dimensions, is escaping into higher dimensions where there is little resistance. 

' This article wishes to point out that almost all known exact cosmological solu- 

tions of deterministic field equations can be understood in terms of a few simple 
ideas. (Perfect fluids with no equation of state or imperfect fluids are meant to 
be excluded by the phrase deterministic field equations.) The very existence of 
these exact solutions and their relationship to each other depends crucially on 
the underlying structure of the field equations themselves, a fact which is rarely 
appreciated. 

These remarks are confined to spatially homogeneous classical solutions of 
gravitational theories in which the purely gravitational variables have Einstein- 
like equations, i.e., do not involve "higher derivatives". For such theories a well 
developed Lagrangian/Hamiltonian formulation of the field equations helps to 
understand the exact solutions in terms of an elegant geometrical picture. This 
approach, initiated by Arnowit, Deser and Misner [1] for Einstein's equations in 
four dimensions, was developed by Misner [2-4] as a powerful way of viewing the 
dynamics of spatially homogeneous cosmological models. It extends to higher 
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dimensional theories in a natural way, including Brans-Dicke and supersymmet- 
ric variations of Einstein's theory coupled to various matter sources. Henneaux 
[5] gives the details for ten-dimensional supergravity, for example. 



2 Spatially Homogeneous Models 

By choosing the time coordinate lines orthogonal to the spatial sections, the 
rf-dimensional spacetime metric can be expressed at least locally in the form 
[6-9] 

('^^g=-N{tfdt(^dt + gab{t)uj''(^w'' (1) 

on the manifold R x (G/H), where are 1-forms on G/H dual to the 

frame {ea}a=i D (let D = d — I) characterized by structure functions C^bc = 

U!'^{[ea,eb]) on G/H. These are constants in the simply transitive case where 
H is the identity subgroup of the isometry group G and {ca} is a left invariant 
frame on G/H = G (i.e., a basis of the Lie algebra of G), but they are functions 
on the left coset space G/H in the nontrivial multiply transitive case, where H 
is a continuous subgroup. In the first case the matrix g = (gab) is an arbitrary 
positive definite symmetric matrix, while in the second case it must satisfy ad- 
ditional linear constraints [8,9]. Additional spatial symmetries in the first case 
lead to similar constraints; such cases might be called trivially multiply tran- 
sitive as opposed to those "nontrivial" models which have no subgroup acting 
simply transitively on the orbits. 

The choice of lapse function N{t) determines the parametrization of the 
family of homogeneous spatial sections, i.e., the time fimction. The proper time 
r, defined by dr = N{t)dt — uj^ corresponds to unit lapse. In expanding models 
with an initial singularity, r is usually chosen so that this singularity occurs at 
r = 0. Defining = dt leads to the spacetime 1-forms oj" dual to the frame 
{ea}a=o,i,...,-D, where eg = d/dt and e± — N^^cq is the unit normal to the 
spatial sections. Clearly the choice of time function should be made to simplify 
the equations of motion for g = (gab)- In order to obtain exact solutions (at least 
up to quadratures) one must decouple the equations in some way, a problem 
which depends critically on the choice of time. 

For the sake of brevity, consider only models which are "diagonalizable" and 
therefore can be assumed to be diagonal, i.e., the matrix g is diagonal 

g = Sd diag(gii , . . . , gnn) (2) 

and hence its natural logarithm is an arbitrary diagonal matrix which can be 
decomposed into its pure trace and tracefree parts 

gD = e'"e'^ , 

ilngD = al + /J^eA , (3) 
TreA = 0, TreAes = £>(£>- , 

where the normalization of the tracefree basis is the most natural one in view 
of considerations described below. Most known exact solutions fall into this 
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diagonal class. The existence of diagonal solutions of the field equations depends 
on the symmetry type, the choice of frame and the form of the field equations 
(including the symmetry of the source, etc.) and is usually connected in some 
way with discrete or continous spacetime symmetry of some kind. 

Even though the present remarks will be confined to such diagonal solutions, 
a similar analysis of the few exact solutions not of this form is even more crucial 
due to their complexity. Moreover, it is absolutely essential to rely on these 
ideas to appreciate the general spatially homogeneous case, which does not 
admit exact solution [10,11]. 

3 Field Equations 

Since we wish to focus on the gravitational properties of higher dimensional 
models, we consider the Einstein equations with all other fields lumped into the 
energy- momentum tensor (including a possible cosmological constant). This 
assumes that the gravitational part of the classical field equations can be repre- 
sented in this way, thus excluding "higher derivative" theories. These equations 

O = M^0^ I^'^^gl'/^^'^^G^IS - kT^p) , l^^^^ = Ng^^ (4) 

(where metric quantities without the leading superscript refer to the spatial 

metric) naturally split into evolution equations = Mab which evolve the spatial 
metric and constraints on the solutions of those equations 

= H = 27V~^M-'-_L = 2M°_L (super-Hamiltonian constraint) , 
= Ha = 2N~^M-^a = 2M°(j (supermomentum constraint) . 

The Einstein equations may also be written in "Ricci form" 

= Pais = ^'^^Raff - kE„0 , = T„0 -{d- 2)-^T\ g^^ , (6) 

a form very convenient for static solutions in discussions of spontaneous com- 
pactification. However, for dynamical solutions a compromise is preferable, 
namely the Ricci evolution equations = Pab, supplemented by the Einstein 
constraints rather than more equations involving second time derivatives. For 
our purposes it is crucial to choose the lapse function to be an explicit function 
of the spatial metric given by 

iVTTG = mo - 1)5'/' ; (7) 

this choice of time will be called Taub time gauge. In fact, the Ricci evolution 
equations are just the variational equations obtained from the usual ADM Ein- 
stein Lagrangian [12] with this explicit choice of the lapse, although for some 
symmetry types one must add a nonpotential force to the system to get the 
correct equations. For diagonal models, these equations are explicitly 

(Ingn)" =-27V|TG(R-fcE) , (8) 
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where R = {R°'b) is the matrix of mixed components of the spatial Ricci tensor 
and E = {E^b). 

In general, the evolution equations may be obtained from the ADM La- 
grangian or its associated Hamiltonian 

-^ADM = T —U , Haum =T + U = NH , 

U = U^ + fjsource ^ V ) 

f/G = N{-g^/'^R) , [/source ^ N{-2g^/'^kT^s_) , 
supplemented by a nonpotential force [6] defined by 

-diN-'U) + Q'^'^dQab = -g^'^G''' - kT'''')dgab (10) 

for some symmetry types, where the exterior derivative is the one on the space 
A4 of metric matrices, on which the DeWitt [13] metric is a Lorentz metric 

2N-'g = 27V-l5V2(3a(y)6 _ gabgcd^dg^^ ^ ^^^^ (^j) 

Allowing the lapse to have a factor depending explicitly on the spatial metric is 
equivalent to conformally rescaling this metric by that factor [4]. In particular, 
the Taub time gauge removes the metric determinant factor and makes {a, (3^} 
natural Lorentz orthonormal coordinates on the flat diagonal submanifold Md 
with respect to the rescaled DeWitt metric (a conventional factor of 2 is missing 
from the kinetic energy) 

2N:^Qg\Mj, = -da (g)da + SAsdp^ O dp^ . (12) 

a is the natural time variable on the space, and (3^ are natural spatial coordi- 
nates. 

If N is considered as an independent variable, one obtains the Einstein evolu- 
tion equations with the super-Hamiltonian constraint as its Lagrange equation. 
Allowing the lapse to depend explicitly on the spatial metric leads to new evo- 
lution equations which differ by multiples of that constraint. (The Hamiltonian 
approach widens this freedom to allow explicit dependence on the gravitational 
momenta.) For diagonal metrics, the useful lapse choices are power law lapses, 
namely products of powers of the diagonal metric components [14,15]. The Taub 
time gauge choice is an example. 

4 Exact Solutions 

All known exact spatially homogeneous solutions of gravitational field equations 
are such that the nongravitational variables decouple from the field equations for 
the gravitational variables in the following sense. The source energy-momentum 
(including a possible cosmological constant) may be represented in terms of the 
metric and constants of the motion, leading to an entirely geometric system. 
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This occurs only when the source has cither discrete or continuous additional 
symmetry. For the diagonal case, the spatial curvature potential is a linear 
combination of individual potential terms which are power law in the diago- 
nal metric components and therefore exponentials whose arguments are linear 
functions of the natural inertial coordinates [6]. To admit exact solutions, the 
source can add at most more such terms to the total potential energy U when re- 
expressed in terms of possible constants of the motion. Independent of this, one 
must have a suitable decoupling of the gravitational modes for exact solutions 
to exist. At most one power of any exponential argument can appear in the to- 
tal potential, i.e., the exponentials with different arguments must be associated 
with linearly independent directions in the space A^d- If this is so, one may 
obtain integrable decoupled equations, the only complication being the causal 
nature of the level hypersurfaces (hyperplancs) of the individual exponential 
potentials. These hyperplanes may be timelike, so that level surfaces move with 
speed ds/da <l'ui ^-^ space, or spacelike or null, in the latter case moving with 
speed ds/da = 1. If different powers of a given exponential argument appear, 
the resulting integrals for decoupled modes become more complicated. 

Consider first the four-dimensional vacuum case (with zero cosmological con- 
stant). The diagonal exact solutions fall into two overlapping families of solu- 
tions which intersect at the well known Kasner family [16]. The first is the family 
found by Taub [17] in his introduction of the spatially homogeneous cosmological 
models. He used the time gauge N = g^^"^, namely the Taub time gauge apart 
from a constant. The second family is due to Joseph [18] and Ellis and MacCal- 
lum [19] and contains a subfamily which can be analytically continued to the 
Kantowski-Sachs models [20]. These latter models are the only nontrivial mul- 
tiply transitive symmetry type for this dimension in the sense that they admit 
no subgroup acting simply transitively on the spatial sections and so cannot be 
represented as a simply transitive case with higher symmetry. The Taub family 
contains one or two nonnuU potentials while the Joseph-EUis-MacCallum family 
contains one null potential in Taub time gauge. (Changing the time gauge may 
change the causality properties.) The Joseph- Ellis-MacCallum solutions have 
been given in the time gauge = (533)^^^ but can easily be converted to the 
Taub time gauge. These two families of exact solutions are very useful examples. 

The Taub family of solutions may be described by 

= = n(i) , C^2 = n(3) , go = e2«e2(/'+e++'5-e-) , ^^^^ 
e+ = diag(l, 1, -2) , e_ = \/3diag(l, -1, 0) , p- = 6{n^'^)Po , 

and their dynamics in Taub time gauge is described by the Hamiltonian (ex- 
pressed in velocity phase space) 

/fTTG = i(-a' + /3+'+r') + 6n(3)V(«-2''^)-2W3)nWe-2('5^-2«) . (14) 
The Lorentz transformation [4] 

{a,P+) = 3-'/\2a-P+,-a + 2p+) (15) 
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leads to a completely decoupled Hamiltonian 



Httg = i(-s2+;^+2 + /j-2) + 6n(3)2e-*^+ -24n(3)n(i)e2^« ,,,, 

(16) 

= -Ha + H-^+ + Hjj- = 

consisting of three 1-dimensional scattering problems with exponential poten- 
tials and constant energies, restricted only by the constraint that the appropri- 
ately signed sum of the individual energies vanish. The solution of an exponen- 
tial scattering problem leads to hyperbolic or trigonometric functions 

= ^{x^+ i^e"^) = E, ^ t = j dx{2E, - iie"^)-'/^ , 
or defining 7 = {2E^y/^ , (17) 



e 2' 



hux _ f m'/'7"' cosh(i7'^i) /i > , E,>Q 
lAtli/27-1 sinhl^-/\iyt\) M < , E^ G ^ . 



(a, /3+./3^) arc incrtial coordinates of the rest frame of the n^^^^ potential, 
which moves with speed d(3'^ /da = | in space. For the Bianchi type IX 
case, characterized by n^^^n^^^ > (so /3~ = 0), the hyperbolic cosine solution 
is relevant, interpolating between the asymptotic free positive and negative ex- 
ponential solutions (for the metric components) at t = ±00; the unit velocities 
in (3^ space of the asymptotic solutions are related by a simple reflection in 
the rest frame of this potential [4]. Letting n'^^^ — > contracts the group to 
Bianchi type 11, eliminating the "tachyonic" n^^-*n*^^) potential (it moves with 
speed d/3~^ /da = 2 in f3^ space) and allows free motion parallel to the n'^^^ po- 
tential (the a and directions). For the Bianchi type Vlll case, characterized 
by n^^^n*^^) < 0, only the positive energy solutions are relevant for the a motion 
(the potential is negative) due to the super-Hamiltonian constraint. When both 
n^^-* and n^'^-' vanish, one obtains the Abclian case of Bianchi type I and all three 
variables have free motion, leading to the Kasner solution when re-expressed in 
proper time gauge, the exponentials in the Taub time converting to powers of 
the proper time. 

The Joseph-EUis-MacCallum family is also described by a Hamiltonian sys- 
tem for a diagonal metric matrix but with structure constants and Taub time 
gauge potential given by 

C^3i=a + q, C^i2 = a-q ,0" = -hq^ , \ = qa'^ = ±{-h)-^'^ , 
Uttg = 12 • 66^^""+^^^ {a^ + gV3) = 12 • 6a2e*(^''+^^) (1 + X^S) . ^^^^ 

The null hyperplane level surfaces of this potential require a different transfor- 
mation of coordinates when a ^ 

-^ = ^0 + f3+ f3° = {l + CysW - P+ + CP- /\/3 

fi+ = f3+ +cf3-/V3 /3+ = p+-cp-/vs-epvs 

f3- = f3- -c{p°+f3+)/v^ r = r+c/37V3, (19) 
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where (J = [1 — 5{a)]X, leading to 

Httg = -i(l + CVmpy ' + 5(^°)' (^+)' + hi^'Y ' + 12 • 6a2e^^°(l + X^S) 
n^3 = -12ae^'^(p+y [note: a;5+ = a/?+ + g/3-/\/3] . 

(20) 

If a 0, then (/?"'")' = (set its constant value equal to zero since it is 
gauge) and the Hamiltonian becomes 

Httg = -(1 + AV3)i?^o + H^- , (21) 

thus decoupling into free motion for 0^ with constant energy (linear behavior) 
and a one-dimensional scattering problem for in a negative exponential po- 
tential with constant energy (the "In sinh" solution) which must be positive to 
satisfy the vacuum super-Hamiltonian constraint of vanishing total energy. 

If a = and therefore ( — 0, then /3+ and /3~ are unchanged by the trans- 
formation while = /S" -|- and it is convenient to introduce another null 
coordinate /3+ = /3+ — /3° so that the kinetic energy becomes 

TTTG = hw'ym- ■ (22) 

Here the null geometry is essential. /3~ and undergo free motion while 

{f3~^)" = 8?7ttg leading to exponential behavior, describing the Bianchi type 
VIq. However, if ^ 0, then the supermomentum constraint requires that 
be a constant (set it equal to zero since it is gauge) while if = one arrives 
again at the Abelian case and all three variables are free. The Bianchi type VIq 
solution in proper time time gauge was given by Ellis and MacCallum [19]. 

The condition a = q ^ (or a = —q) describes Bianchi type III = VI_i. The 
vacuum supermomentum constraint then imposes the condition (.gn/gss)' = 
associated with local rotational symmetry, the spatial geometry being that of a 
family of geodesically parallel 2-dimensional surfaces of constant negative cur- 
vature. Allowing to take negative values in the Einstein equations simply 
changes the sign of the potential and describes the Kantowski-Sachs models 
which instead have constant positive curvature on the geodesically parallel fam- 
ily of 2-dimensional surfaces, changing the solutions of the one-dimensional scat- 
tering problem from hyperbolic to trigonometric functions. 

All of these solutions describe the effect of zero, one or two independent 
exponential potentials on the flat diagonal configuration space A4b- They can 
be extended to certain nonvacuum nonzero cosmological constant cases in four 
and higher dimensions in a trivial way by a mechanism called "variation of pa- 
rameters" [15]. The constraints impose certain values on the parameters which 
appear in the evolution equations. Allowing the parameters to vary from their 
constrained values is equivalent to introducing a nonzero energy-momentum 
tensor which may be interpreted as that of a stiff perfect fluid moving orthog- 
onal to the spatial sections if the supermomentum is zero, or equivalently as a 
massless scalar field. Nonzero values of the supermomentum in the Joseph-Ellis- 
MacCallum case lead to tilted fluid flow and can also be interpreted in terms 
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of changing the spacctimc symmetry from isometries to homothetic transforma- 
tions, leading to spatially self-similar solutions. 

The massless scalar field can be obtained by dimensional reduction from a 
5-dimcnsional vacuum Kaluza-Klcin theory, together with a Wcyl transforma- 
tion, namely a conformal transformation of the four-dimensional metric when 
imbedded in five dimensions [21]. This allows one to relate the above solu- 
tions to solutions of the Brans-Dickc theory [21-23]. One may also extend the 
imbedding to vacuum Kaluza-Klein theories with any number of additional flat 
dimensions, thus obtaining exact solutions of those theories [7,24]. The Weyl 
transformation is related to a redefinition of the variable a when adding extra 
dimensions. 

Nonvacuum Kaluza-Klein theories also admit such an extension of the four- 
dimensional vacuum solutions by variation of parameters. Both ten and eleven- 
dimensional supergravity when restricted to the Bose sector merely add ad- 
ditional exponential potentials to the spatial curvature potentials when one 
assumes a Frcund-Rubin [25] or Freund-Rubin-Englert [26] ansatz, since the 
differential form fields can then be expressed in terms of the metric and con- 
stants of the motion. The scalar field in the ten-dimensional case is related by 
dimensional reduction to an extra flat dimension [27]. One thus sees the same 
familiar exact solutions reappearing [28-31]. 

These solutions may also be deformed by a second kind of variation of param- 
eters in which one changes the parameters which appear in the field equations 
themselves in such a way that the equations still remain integrable. In the 
Abelian case, the addition of a cosmological constant or a perfect fluid with 
equation of state p = (j — l)p whose flow is orthogonal to the spatial sections 
adds a single a-dependent exponential to the Hamiltonian, a timelike potential 
like one of the two in the semisimple Taub case in four dimensions. One can add 
such a potential to the nonsemisimple Taub case or the Joseph-EUis-MacCallum 
case and still obtain the same kinds of solutions. On the other hand, adding a 
locally rotationally symmetric electromagnetic field changes the nature of the 
solutions, but only deforms them in a simple way, leading to more very similar 
exact solutions. One may in fact add many fields simultaneously, leading to, for 
example, solutions in four dimensions with a stiff perfect fluid, a scalar fleld and 
an electromagnetic field, the scalar field allowing the optional reinterpretation 
in terms of solutions of the Brans-Dicke theory. 

Consider adding a locally rotational electromagnetic fleld to the Taub family 
of models in four dimensions [32]. This simply adds another /?"*" potential to the 
Hamiltonian but with a different power, leading to a one-dimensional scattering 

problem with two different exponential potentials. The explicit integral in (17), 

_ 1 

when rewritten in terms of the variable u = /i^^"e s''^^ becomes proportional 
to the standard integral ( / du{c^v?' — 1)~^/^) for the inverse hyperbolic cosine 
when fj, and are positive as occurs in this case. The additional potential 
merely adds a linear term to the radical which can be eliminated by completing 
the square. The new solution for u is then a constant plus another constant 
times the old solution, a simple variation of parameters [15,33]. 
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The Taub time gauge is clearly very powerful, but other power law lapse 
time gauges are also useful. Changing the lapse leads to different decoupling 
possibilities for different kinds of driving terms. For example, adding a cosmo- 
logical constant to the semisimplc Taub case (Bianchi types VIII and IX) leads 
to three exponential potentials along distinct directions in the two-dimensional 
Taublike subspace of Mb associated with local rotational symmetry. In addition 
to the orthogonal pair of timelike and spacelike potentials, a cosmological con- 
stant term adds a null potential, causing a coupling of the equations. However, 
in Misner's [15,34] time gauge N ~ (533)"^^^, one degree of freedom remains 
unchanged but enters the equation of motion for the second in such a simple 
way that the new solution merely adds on a term to the old one [35]. There 
are many other instances in various time gauges where similar ideas apply (in- 
cluding the .Joscph-EUis-MacCallum family in the original time gauge and even 
the familiar Friedmann models in conformal time gauge iV ~ e", proper time 
time gauge and other time gauges, as well as higher dimensional generalizations 
of these solutions). However, no systematic analysis of lapse choices has ever 
been undertaken as a way of simplifying the equations of motion, other than 
the limited discussions of references [7] and [15]. 

It is incredible how many pages of journal space have regurgitated the same 
old solutions in various disguised forms without leading to some recognition of 
the connections. Too many people interested in exact cosmological solutions 
suffer to some degree from a sort of tunnel vision which prevents them from 
taking the time to appreciate the ADM perspective for the very branch of grav- 
itation for which it is the most powerful. The present article shows that the 
rewards for doing this are considerable. 
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